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ABSTRACT 

We present predictions for the line-of-sight velocity dispersion profiles of dwarf 
spheroidal galaxies and compare them to observations in the case of the Fornax dwarf. 
The predictions are made in the framework of standard dynamical theory of spherical 
systems with different velocity distributions. The stars are assumed to be distributed 
according to Sersic laws with parameters fitted to observations. We compare predic- 
tions obtained assuming the presence of dark matter haloes (with density profiles 
adopted from TV-body simulations) versus those resulting from Modified Newtonian 
Dynamics (MOND). If the anisotropy of velocity distribution is treated as a free pa- 
rameter observational data for Fornax are reproduced equally well by models with 
dark matter and with MOND. If stellar mass-to-light ratio of IMq/Lq is assumed, 
the required mass of the dark halo is 1.5 x IO^Mq, two orders of magnitude bigger 
than the mass in stars. The derived MOND acceleration scale is gq — 2.1 x 10^^ cm/s^. 
In both cases a certain amount of tangential anisotropy in the velocity distribution is 
needed to reproduce the shape of the velocity dispersion profile in Fornax. 
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1 INTRODUCTION 

Dwarf spheroidal (dSph) galaxies provide a unique testing 
field for the theories of structure formation. Due to their 
large velocity dispersions they are believed to be dominated 
by dark matter (hereafter DM, for a review see Mateo 1997). 
Because of their small masses they also lie in the regime of 
small accelerations, where, according to an alternative the- 
ory of Milgrom (1983a), the Modified Newtonian Dynamics 
(MOND) may operate causing a similar increase of velocity 
dispersions. 

Until recently the measurements of velocity dispersion 
in dwarfs were restricted to a single value, the central veloc- 
ity dispersion, on which all of the dynamical modelling was 
based. In order to estimate the mass-to-light ratio in such 
case a simplifying assumption of isotropic velocity distribu- 
tion was usually made. Using this approach many authors 
(see Mateo 1998 and references therein) estimated mass-to- 
light ratios in dwarfs to be much larger than stellar val- 
ues and interpreted this as a signature of the presence of 
DM. However, Milgrom (1995) analysed the data concern- 
ing seven well known dSph galaxies and concluded that when 
MOND is applied their mass-to-light ratios are consistent 
with stellar values if all the observational errors are prop- 
erly taken into account. 

Given the velocity dispersion data as a function of dis- 



tance, which are becoming available now, we are able to 
extract much more information from the observations. Our 
conclusions do not have to rely on the assumed velocity 
anisotropy model, but the best-fitting model can actually 
be found. On the theoretical side some progress has also 
been achieved. The DM distribution in haloes can now be 
reliably modelled using the results of A'^-body simulations 
(e.g. Navarro, Frenk & White 1997, hereafter NFW). 

The purpose of this work is to explore the two hypothe- 
ses, based on DM and on MOND, by calculating the pre- 
dicted velocity dispersion profiles and comparing them to 
available data. In the case of DM our analysis is somewhat 
similar to an earlier study of Pryor & Kormendy (1990) who 
also considered two-component models involving stars and 
DM. In the case of MOND, our approach here is related 
to that of Milgrom (1984) who applied the Jeans equa- 
tion to study density profiles of equilibrium systems with 
isotropic and isothermal velocity distributions. Here we as- 
sume instead the observed density distribution of stars to 
obtain the velocity dispersions in MOND for different veloc- 
ity anisotropy models. 

2 VELOCITY DISPERSION PROFILE 

The radial velocity dispersion (Jr{r) of stars can be obtained 
by solving the Jeans equation (Binney & Tremaine 1987) 
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dr r 



— ua^ = vg, 



(1) 



where v{r) is a 3D density of stars, g{r) is the gravi- 
tational acceleration associated with the potential $(r), 
g — — d<l>/dr, and /3 = 1 — ag{r)/a^{r) is a measure of the 
anisotropy in the velocity distribution. 

In the absence of direct measurements of velocity 
anisotropy in dSph galaxies we have to consider different 
models for /3. Our choice here is guided by the results of 
Ai'-body simulations of DM haloes and observations of el- 
liptical galaxies. For DM haloes Thomas et al. (1998) find 
that, in a variety of cosmological models, the ratio agja-c is 
not far from unity and decreases slowly with distance from 
the centre to reach ~ 0.8 at the virial radius. Observations 
of elliptical galaxies are also consistent with the orbits being 
isotropic near the centre and somewhat radially anisotropic 
farther away, although cases with tangential anisotropy are 
also observed (Gerhard et al. 2001). 

Our fiducial anisotropy model will therefore be that of 
isotropic orbits: ae(r) = o^ir) and /3 = 0. A more realistic 
approximation is provided by a model proposed by Osipkov 
(1979) and Merritt (1985) with /3 dependent on distance 
from the centre of the object 



/3oM = 



+ r-i 



(2) 



where r^ is the anisotropy radius determining the transition 
from isotropic orbits inside to radial orbits outside. This 
model covers a wide range of possibilities from isotropic or- 
bits (ra — » oo) to radial orbits (ra — > 0). Another possibility 
is that of tangential orbits with the fiducial case of circular 
orbits when /? -^ —oo. To cover intermediate cases between 
the isotropic and circular orbits we will consider a simple 
model of /3 =const< 0. 

The solution of the Jeans equation (|l|) with the bound- 
ary condition cTr — > at r — » cxa for /3om anisotropy is 



ua^{P = /3om) 



1 



(r +r^)ug dr. 



while for /3=const we find 

ua,. {/3 — const) — —r~ j r vg dr. 



(3) 



(4) 



Prom the observational point of view, an interesting, 
measurable quantity is the line-of-sight velocity dispersion 
obtained from the 3D velocity dispersion by integrating 
along the line of sight (Binney & Mamon 1982) 



o-L(-R) = 



I{R) 



R^ 



dr 



(5) 



where I{R) is the surface brightness. For circular orbits, 
(J, = 0, and one has 



^los{R) = 



I{R) 



(?) 



uV^r 



^r 



i?2 



dr 



(6) 



where V = {rgY'^ is the circular velocity associated with 
the gravitational acceleration g{r)- 

Introducing results (|3|) or 1^ into equation (|5|) and in- 
verting the order of integration the calculations of aios can 
be reduced to one- dimensional numerical integration of a 
formula involving elementary functions in the case of /3om 



(see Mamon & Lokas, in preparation) and special functions 
for arbitrary (3 =const. 

We will now discuss difi^erent factors entering in the cal- 
culations of aios, namely the density distribution of stars v, 
the mass distributions contributed by stars and DM to grav- 
itational acceleration g and how this acceleration is changed 
in MONO. 

2.1 Stars 

The distribution of stars is well known in a number of dSph 
galaxies from surface brightness measurements and is usu- 
ally well fitted by the exponential law or its generalization 
i.e. the Sersic profile (Sersic 1968, see also Ciotti 1991) 



7(i?) = /oexp[-(i?/i?s 



(7) 



where 7o is the central surface brightness and Rs is the char- 
acteristic projected radius of the Sersic profile. The Sersic 
profile with parameter m varying in the range 1 < m < 10 
has been found (Caon, Capaccioli & D'Onofrio 1993) to fit 
surface brightness profiles of elliptical galaxies in large mass 
range better than the standard de Vaucouleurs law (corre- 
sponding to m = A). For dSph systems the fits are usually 
performed with m = 1, although is some cases other values 
of m are found to provide better fits (e.g. Caldwell 1999) . 

The 3D luminosity density v{r) is obtained from I{R) 
by deprojection 



Kr) 



Al AR 



TT y,, AR ^R2 _ r2 



(8) 



In the case of m — 1 we get ^{r, m = 1) = 
IoKo{r/R3)/{nR3), where Ko{x) is the modified Bessel 
function of the second kind. For other values of m in the 
range 1/2 < m < 10 an excellent approximation for z/(r) is 
provided by (Lima Neto, Gerbal & Marquez 1999) 



^{r) 



1^0 



/or(2m) 



/ J, \ l/m 

[r^j 



(9) 



2Rsr[{S - p)m] 
= 1.0 - 0.6097/m -f- 0.05463/m^ . 
The mass distribution of stars following from (H) is 

7[(3-p)m,(r/i?s)^/'"] 



M,(r) = rL 



(10) 



r[(3-p)m] 

where T=const is the mass-to-light ratio for stars, Ltot is the 
total luminosity of the galaxy and "f{a, x) — J e~*i°'~^di 
is the incomplete gamma function. 

2.2 Dark matter 

We assume that the density distribution of DM is given by 
the so-called universal profile (NFW, see also Lokas & Ma- 
mon 2001) 

Pcrit,o (r/rs) (1 -f r/rs)2 ^ ' 

with a single fitting parameter 5char, the characteristic den- 
sity. The scale radius rg is defined by r^ = r^/c where r^ is 
the virial radius (the distance from the centre of the halo 
within which the mean density is 200 times the present crit- 
ical density, pcrit,o) and c is the concentration parameter 
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related to the characteristic density by Schm = 200c^(?(c)/3 
with g{c) = l/[ln(l + c) — c/(l + c)]. The value of c is known 
to depend on mass and will be extrapolated here to small 
masses characteristic of dwarf galaxies assuming the follow- 
ing formula fitted to A''-body results for ACDM cosmology 
of Jing & Suto (2000) 



c= 10.23 



1012 Mq 



(12) 



where Mv = 8007rrvpcrit,o/3 is the mass within virial radius 
and h is the Hubble constant in units of 100 km s~i Mpc~i. 
In the following we will use h = 0.7. 

The mass distribution associated with the NFW profile 
©is 



MNFw(r) = g{c)M^ 



ln(l + cs) 



cs 



1 + cs 



(13) 



where s — r/r^, and for h = 0.7 the virial radius scales with 
the virial mass so that r^ = 206[Afv/(10^^AfQ)]i''^ kpc. 

2.3 MOND 

In the case of MOND the gravitational acceleration is 
changed so that (Milgrom 1984) 



5m 



9N 



K9N/ao)' 



(14) 



where qm and qn are the modified and Newtonian gravi- 
tational accelerations, respectively, ao is the characteristic 
acceleration scale of MOND estimated to be of the order of 
10~* cm/s'^. The function fJ.{x) has to reproduce Newtonian 
dynamics at large accelerations, i.e. i-i{x) -^ 1 for x — > oo, 
and is supposed to be i-i{x) ~ x for a; -^ 0. A good ap- 
proximation is then provided by e.g. fJ.{x) = x{l + x'^)'^'"^ 
(Milgrom 1983b). Using this approximation we find from 
equation (fl4|) 



flM 



ffN+ffN(4ao-h5N 



2N1/2 



1/2 



(15) 



3 RESULTS FOR THE FORNAX DWARF 

The structural parameters of the Fornax dwarf needed in the 
following computation were taken from Irwin & Hatzidim- 
itriou (1995, hereafter IH). They fitted the surface brightness 
profile of Fornax with the exponential distribution (m = 1 in 
equation |7|) finding the Sersic radius (re in their notation) 
Rs ~ 9.9 arcmin = 0.35 kpc assuming the distance modu- 
lus m — M = 20.4. The brightness of Fornax in V-band is 
My = —13.0 mag so the total luminosity needed in equation 



(0) is Lt, 



— 1.4x 10 Lq. We also adopt the mass-to-light 
ratio for stars in this band to be Ty ~ 1Mq/Lq (Mateo et 
al. 1991). 

Figure 111 shows the predicted line-of-sight velocity dis- 
persion profiles for the Fornax dwarf together with the data 
from Mateo (1997). The two panels present results for differ- 
ent velocity anisotropy models. The upper panel gives the 
results in the case of isotropic orbits obtained from equa- 
tions (|3|) and (g) in the limit of large ra. The lower panel is 
for circular orbits showing results from (pi). The three types 
of curves in each panel correspond to the three possibilities 
of the velocity dispersion profile being generated 1) by stars 
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Figure 1. Upper panel: predicted line-of-sight velocity dispersion 
profiles for the Fornax dwarf in the case of isotropic orbits. The 
dotted line shows the prediction obtained by solving Jeans equa- 
tion for stars only. The three solid lines correspond to stars -|- DM 
with the dark halo masses of 10®, 10® and lO^" A/q from bottom 
to top. The three dashed lines show results for stars in MOND 
with ao = 1, 3, 6 xlO~* cm/s^ from bottom to top. The data 
for Fornax are from Mateo (1997). Lower panel: same as upper 
panel, but for circular orbits. 

moving under Newtonian gravitational acceleration gener- 
ated by stars only (dotted curves), 2) by stars moving un- 
der Newtonian gravitational acceleration generated by stars 
and DM (solid lines), 3) by stars moving under modified 
(MONDian) acceleration generated by stars (dashed lines). 
For all cases the 3D distribution of stars v(r) is given by 
(H). The gravitational acceleration for the Newtonian cases 
1) and 2) is g = pN = -GM{r)/r'^ with M{r) = Mt{r) for 
1) and M{r) = Af*(r) -I- MNFw(r}for 2), where Af*(r) and 



r) fo: 

m 



AfNFw(r) are given by equations ( llCJ ) and (|l^), respectively. 
In the case of MOND we have g = gm from equation ( |l5[ ) 
with (7n = —GMf{r)/r'^ . The three solid curves in each panel 
correspond from bottom to top to the virial masses of the 
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Figure 2. Best-fitting models for Fornax based on stars+DM 
(solid line) and MOND (dashed line). The model involving DM 
has /3 = -1.4 and Mv = 1.5 X IO^Mq. The model with MOND 
requires f) = —1 and ag = 2.1 X 10~* cm/s^. 

dark halo Mv = 10*, 10^ and 10^° Mq. The three dashed 
lines show from bottom to top results obtained with different 
characteristic MOND acceleration scales ao = 1, 3, 6 xlO~* 
cm/s^ . 

Figure hi clearly shows that the stars alone cannot re- 
produce the observed velocity dispersion profile of Fornax, 
unless the stellar mass-to-light ratio is much larger than as- 
sumed. In the following we will farther analyse the cases 
2) and 3) trying to find the best-fitting models. Compar- 
ing the predictions shown in Figure ry with the data we find 
that a certain amount of tangential anisotropy is needed 
to reproduce the shape of the observed velocity dispersion 
profile. Indeed, trying more radial anisotropy than in the 
case of isotropic orbits by using_small values of ra (compa- 
rable to Rs or r\) in dd) and (H) makes the curves in the 
upper panel of Figure |l[decrease even more steeply with dis- 
tance, contrary to the trend observed in the data. Therefore 
in what follows we adopt the velocity anisotropy model of 
/3 =const< and treat /? as a free parameter. 

To find the best models we performed a two-parameter 
fitting for each case. The fitting parameters were (Mv,/3) 
for stars-l-DM case and (ao,/?) for MOND. The resuhs of 
minimizing x^ weighted by the variances of the data points 
are shown in Figure 0. It turns out that the DM model fits 
the data best for /3 = -1.4 and Mv = 1.5 x lO'' Mq. This 
virial mass is more than two orders of magnitude larger than 
the stellar mass (1.4 x IO'^Mq for Tv = 1Mq/Lq). In the 
case of MOND the best-fitting parameters are /3 = — 1 and 
ao = 2.1 X 10~* cm/s'^. This value of ao is consistent with 
other estimates (Milgrom 1983b). 

Figure S shows the combined mass-to-light ratio in V- 
band for Fornax calculated from 
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(16) 
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where Mnfw(»") and Af, (r) are given by equations (|l3|) and 



Figure 3. The dependence of the mass-to-light ratio of Fornax 
on distance. The mass comes from stars and DM with the best 
fitting Afv = 1.5 X IO'^Mq. 

(iLy) respectively, and Ly{r) — M, (r)/Tv is the luminosity 
distribution. The shape of the curve refiects the fact that 
DM dominates both at small and large distances. The di- 
vergence of M/L at the centre of the galaxy is due to the 
fact that the luminosity density (Mj increases towards the 
centre more slowly than the DM cusp of r^^ in (|ll|). The 
growth of M/L at large distances is due to the behaviour 
of AfNPw(r) which diverges logarithmically while the mass 
in stars is finite. The minimum of M/Lv ~ 11 appears at 
a scale of the order of the Sersic radius and is consistent 
with the values of the central mass-to-light ratio estimated 
by Mateo et al. (1991). 

4 DISCUSSION 

We have found that the predictions based on MOND with 
ao = 2.1 X 10~* cm/s'^ fit the Fornax data equally well as 
models based on DM with virial mass of A/v ~ 1.5 x 10^ Mq. 
In both cases a certain amount of tangential anisotropy is 
needed to reproduce the shape of the observed velocity dis- 
persion profile. It turned out that the MOND predictions 
in the case of Fornax dwarf need less tangential anisotropy 
than DM predictions to reproduce the shape of the velocity 
dispersion profile. It remains to be checked how the DM pre- 
dictions depend on the particular shape of DM distribution. 
The particular fitted values of ao and A/v are of course 
sensitive to the observational parameters assumed. The mea- 
sured luminosity, distance and structural parameters of the 
dwarfs are all known with limited accuracy, but the de- 
tailed analysis of the errors in them is beyond the scope 
of this paper. In the case of MOND the most doubtful in- 
put parameter is the stellar mass-to-light ratio in V-band. 
We have assumed here Ty ~ 1Mq/Lq characteristic for 
a relatively young stellar population of Fornax, however 
this value may well be different by a factor of two (Ma- 
teo et al. 1991). To estimate its effect we performed a simi- 
lar MOND fitting as described in the previous section with 
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Figure 4. Best-fitting models for Draco (left panel) and Ursa 
Minor (right panel) based on stars+DM (solid lines) and MOND 
(dashed lines). The data are from Armandroff et al. (1997) 



Tv = 0.5 and 2 Mq/Lq, obtaining respectively ao — 4.3 
and 1.1 X 10~* cm/s^ with the same best-fitting /9. If 
Tv = 0.5Mq/Lq is indeed the lower limit of mass-to-light 
ratio, then ao = 4.3 x 10~* cm/s^ is the upper limit for ao. 

In the case of DM the stellar mass-to-light is not likely 
to influence much the results since the input from stars 
to the velocity dispersion is small. Here the most doubt- 
ful parameter is the concentration of DM haloes, which we 
adopted to follow the same dependence on mass as larger 
haloes actually studied in A''-body simulations. Given, how- 
ever, that this dependence is weak in the ACDM Universe, 
the possible error associated with this factor is rather small. 

One more reason to worry are the tidal interactions that 
may affect the dynamics of some dwarfs and our interpre- 
tations of their velocity dispersions. Piatek & Pryor (1995) 
have performed numerical simulations of such effects and 
concluded that in Newtonian dynamics they should not af- 
fect much the inferred mass-to-light ratios (but see Kroupa 
1997). However, Brada & Milgrom (2000) found that in 
MOND, due to its non-linearity, a galaxy system can be 
more affected by an external field. 

One may ask if similar results hold for other dSph galax- 
ies as well. Of those studied by IH the velocity dispersion 
profiles with a few data points are available only for Sex- 
tans, Ursa Minor and Draco dwarfs (Hargreaves et al. 1994a, 
1994b, 1996; Armandroff, Olszewski & Pryor 1995; Arman- 
droff, Pryor & Olszewski 1997). Applying the same analysis 
as we did for Fornax to the 5 data points for Draco and Ursa 
Minor from Armandroff et al. (1997) (with Ty = IMq/Lq) 
we find that again the MOND and DM models can fit the 
data equally well. For Draco the best-fitting models have 
tangential anisotropy, we get /? = —0.9, Mv = 4.2 x 10® Mq 
(2 X 10* the mass in stars) for DM and /3 — —1.5, ao = 
50 X 10"** cm/s^ for MOND. For Ursa Minor the best fits re- 
quire a certain amount of radial anisotropy which we model 
with (|). We get r^ = O.OUrv = OMRs = 0.16 kpc, 
Mv = 1.9 X 10* Mq (10^ the mass in stars) for DM and 
ra = 0.012rv = 0.76i?s = 0.14 kpc, ao = 12 x 10"* cm/s^ 
for MOND. The predictions of the best-fitting models to- 
gether with the data are shown in Figure |j. 

The case of Ursa Minor is doubtful since it has been 
recently observed to possess tidal tails (Martmez-Delgado 
et al. 2001) and therefore may not be in virial equilibrium. 
The high value of ao estimated for Draco can, however, pose 



a serious problem for MOND (we agree on this point with 
Gerhard & Spergel 1992) since ao is supposed to be a univer- 
sal constant. Bringing the value down by an order of mag- 
nitude would require very high stellar mass-to-light ratio 
Tv ~ 1OMq/I/0 or we would have to believe that the mea- 
surements seriously overestimate the velocity dispersion. 

ACKNOWLEDGEMENTS 

I wish to thank A. Graham, P. Kroupa, G. Mamon, C. Pryor, 
P. Salucci and the referee, M. Mateo, for useful comments 
and suggestions which helped to improve the paper. I am 
also grateful to Y. P. Jing for providing the results for the 
concentration of DM haloes in numerical form. This research 
was partially supported by the Polish State Committee for 
Scientific Research grant No. 2P03D02319. 

REFERENCES 

Armandroff T. E., Olszewski E. W., Pryor C, 1995, AJ, 110, 2131 

Armandroff T. E., Pryor C, Olszewski E. W., 1997, in Sofuc Y., 
ed., Proc. lAU Symp. 184, The Central Regions of the Galaxy 
and Galaxies. Kluwer, Dordrecht, p. 35 

Binney J., Mamon G. A., 1982, MNRAS, 200, 361 

Binney J., Tremaine S., 1987, Galactic Dynamics. Princeton Univ. 
Press, Princeton, chap. 4.4.3 

Brada R., Milgrom M., 2000, ApJ, 541, 556 

Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 1013 

Caldwell N., 1999, AJ, 118, 1230 

Ciotti L., 1991, A&A, 249, 99 

Gerhard O., Kronawitter A., SagUa R. P., Bender R., 2001, AJ, 
121, 1936 

Gerhard O. E., Spergel D. N., 1992, ApJ, 397, 38 

Hargreaves J. C, Gilmore G., Irwin M. J., Carter D., 1994a, MN- 
RAS, 269, 957 

Hargreaves J. C, Gilmore G., Irwin M. J., Carter D., 1994b, 
MNRAS, 271, 693 

Hargreaves J. C, Gilmore G., Irwin M. J., Carter D., 1996, MN- 
RAS, 282, 305 

Irwin M., Hatzidimitriou D., 1995, MNRAS, 277, 1354 

Jing Y. P., Suto Y., 2000, ApJ, 529, L69 

Kroupa P., 1997, NewA, 2, 139 

Lima Neto G. B., Gerbal D., Marquez I., 1999, MNRAS, 309, 481 

Lokas E. L., Mamon G. A., 2001, MNRAS, 321, 155 

Martinez-Delgado D, Alonso-Garcia J., Aparicio A., Gomez- 
Flechoso M. A., 2001, ApJ, 549, L63 

Mateo M., 1997, in Arnaboldi M. et al., eds, ASP Conf. Ser. Vol. 
116, The Nature of Elliptical Galaxies. Astron. Soc. Pac, San 
Francisco, p. 259 

Mateo M., 1998, in Richtler T., Braun J. M., eds, Proc. 
Bonn/Bochum-Graduiertenkolleg Workshop, The Magellanic 
Clouds and Other Dwarf Galaxies. Shaker Verlag, Aachen, p. 
53 

Mateo M., Olszewski E. W., Welch D. L., Fischer P., Kunkel W., 
1991, AJ, 102, 914 

Merritt D., 1985, AJ, 90, 1027 

Milgrom M., 1983a, ApJ, 270, 365 

Milgrom M., 1983b, ApJ, 270, 371 

Milgrom M., 1984, ApJ, 287, 571 

Milgrom M., 1995, ApJ, 455, 439 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 

Osipkov L. P., 1979, PAZh, 5, 77 

Piatek S., Pryor C, 1995, AJ, 109, 1071 

Pryor C, Kormendy J., 1990, AJ, 100, 127 

Sersic J. L., 1968, Atlas de Galaxies Australes, Observatorio As- 
tronomico, Cordoba 

Thomas P. A. et al., 1998, MNRAS, 296, 1061 



© 0000 RAS, MNRAS 000, 000-000 



